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Abstract. We obtain a global weighted LP estimate for the gradient of the 
weak solutions to divergence form elliptic equations with measurable coeffi- 
cients in a nonsmooth bounded domain. The coefficients are assumed to be 
merely measurable in one variable and to have small BMO semi-norms in 
the remaining variables, while the boundary of the domain is supposed to be 
Reifenberg fiat, which goes beyond the category of domains with Lipschitz 
continuous boundaries. As consequence of the main result, we derive global 
gradient estimate for the weak solution in the framework of the Morrey spaces 
which implies global Holder continuity of the solution. 



1. Introduction 

This work is concerned with weighted L^-regularity of the gradient of weak 
solutions to the Dirichlet problems regarding elhptic equations with possibly mea- 
surable coefficients in nonsmooth domains. The problems in mind are related to 
some important variational problems arising in the mechanics of membranes and 
films of simple nonhomogeneous materials which form a linear laminated medium. 
In particular, a highly twinned elastic or ferroelectric crystal is a situation where a 
laminate appears. The equilibrium equations of such linear laminates usually have 
merely bounded measurable coefficients, see [3l [TTl [T8[ IT9l [24] . 

In this paper we consider a nonhomogeneous elliptic equation in divergence form 
with bounded measurable coefficients in a very nonsmooth domain beyond the 
class of domains with Lipschitz continuous boundaries. Precisely, we deal with the 
Dirichlet problem 

I A {a^'ix)D,u) - D,f{x) in 

I u — on dfl. 

Throughout the paper the standard summation notation is employed for 1 < 
i,j < n with n > 2. Here J7 is a bounded open domain in with boundary dfl 
and F{x) — (/^(x), . . . , /"(x)) is a vector valued function belonging to a suitable 
weighted space. The matrix of the coefficients A{x) = {a*-' (a:)} : M" — >■ K"^ is 
assumed to be uniformly bounded and uniformly elliptic. Namely, we suppose that 
there exist positive constants L and v such that 

(1-2) ||A||^^(^„^„„.) < L, 
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and 

(1.3) a^H^)^^^J > ^ICP 

for all vectors ^ € M" and for almost every x G M". 

It is well known that under these basic assumptions imposed on A and fi, the 
Dirichlet problem (jl.ip has a unique weak solution if |Fp G L-'^(J7). That is, the zero 
extension -u of u belongs to Hq{W-) and satisfies the corresponding weak integral 
formulation 

(1.4) / a'WjuDiip dx= fD.if dx 
Jn Jn 

for all ip e i?Q(f2). Moreover, the standard i^-estimate 

(1.5) [ \Du\^ dx<c [ \F\^ dx 

Jn Jn 

holds true for the gradient Du of the weak solution with a positive constant c 
depending on n, v, L and the Lebesgue measure of the domain £7. 

A natural extension of the L^-estimate (jl.Sp is the following LP-estimate 

(1.6) / \DuY' dx<c [ \F\P dx 

Jn Jn 

with exponent 1 < p < oo for some c = c(n, i^, L,p, \ > 0. Needless to say, for 
this estimate to be valid, the basic structure requirements (|1.2|) and (ll.3|l on A 
and the boundedness of fl are generally not enough. Some additional regularity 
condition on A and some finer geometric assumption on dfl must be imposed. A 
classical problem is to find such a minimal conditions on A and dQ under which the 
estimate (|1.6p is true for all p in the range (1, oo). This is the so-called optimal W^'^- 
regularity (or equivalently, maximal regularity) problem regarding (jl.ip . As far as 
such minimal conditions on A and for the W^-'^'P-regularity, we are dealing here 
with coefficients matrix A having entries of small bounded mean oscillation (BMO) 
with respect to some of the variables and Reifenberg fiat domains $7, respectively. 
We refer to [201 ES [Ml HI] for the definitions and the basic properties of the BMO 
space and the Reifenberg flat domains (see also [H 13 |33l |34l l35]). 

There have been many research activities on the W^^'^'-regularity problem, cf. [TJ 
[21 HI El ini HSl Uni 133 133 for instance, most of these considering principal coefficients 
of the elliptic operator belonging to the spaces of functions with vanishing mean 
oscillation (VMO) or small BMO. However, little is known about that question in 
the case of only measurable coefficients a*^ . There is a classical by now example 
due to Meyers [26j , showing that if these are merely measurable with respect to two 
independent variables, then the Vl^^'^'-regularity of (jl.ip fails in general. Indeed, 
the Meyers counterexample is easily extendable to the n-dimensional case (n > 2) 
of equations with coefficients that are only measurable with respect to two of the 
variables. 

According to the recent works in [6l one can allow the coefficients a*^ to 
be merely measurable in one of the variables for the estimate p.6p to be true, 
while the boundary of the domain 17 belongs to the class of Reifenberg flat domains 
which goes beyond the category of sets with Lipschitz continuous boundaries. More 
precisely, the global LP-estimate holds true for all p g (1, oo) if for each point and 
for each scale the coefficients are only measurable in one variable and are averaged 
in the sense of small BMO with respect to the remaining n — 1 variables, while the 
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boundary dfl can be trapped into two hyperplanes depending on the scale chosen. 
Let us emphasize on the fact that there is no any regularity assumption with respect 
to one of the variables and the boundary of the domain, being Reifenberg flat, can 
have rough enough fractal structure (see the excellent survey by Toro ^41j). 

Very recently, Dong, Kim and Krylov in [l6l[T7l[22] obtained a global L^-estimate 
under similar assumptions as these in [6l |8]. The approach used there is based 
on making use of a priori pointwise estimates and working on the sharp maximal 
function of the gradient of solutions. This interesting technique is first developed by 
Krylov in [3T^ and gives a unified approach for both divergence and nondivergence 
form equations. 

The present article is a natural outgrowth of [8] and deals with weighted W^ ^- 
theory for the Dirichlet problem (jl.ip . In particular, we derive an extended version 
of the L^'-estimate (|1.6|) in the settings of the weighted Lebesgue spaces 
generalizing this way the M^^'P-regularity theory recently elaborated in [B1IH]. More 
precisely, we prove that under the same regularity assumptions on A and as these 
in [6j|8], the following global weighted Vl^^'^'-regularity 

(1.7) \F\^eLi{n) =^ \Du\''&LUn) 

holds for any p e (2, oo) and with a weight w belonging to the Muckenhoupt class 
Ae (we refer the reader to Section [2] for the precise definitions and notations). Let 
us point out that the W^'P-regularity derived in [HI H] is a special case of (II. 7p when 
w{x) = 1, which makes (|1.7|) a natural extension of the M^^'P-theory. 

Throughout the paper, the case in mind is when 2 < p < oo and the prescribed 
assumption on the free term in is 

\Ff e Ll{n), w £ Ae, 2<p<oo. 

Under this condition, needless to say, the problem (II. ip has a unique weak solution 
u G H^(fl) because Li]{n) C L^{fl). Assuming that at each point and for each scale 
the coefficients of (jl.ip are only measurable in one variable and are averaged in the 
sense of small BMO with respect to the remaining n — 1 variables, and considering 
domains 51 with Reifenberg fiat boundaries, we will prove that 

\Du\^ G Li{n) 
with the corresponding gradient estimate 

(1.8) w\Dufw^,^^^^<4m^i^^^ 

for every p G (2,oo), where the constant c is independent of u and F. 

Similar weighted results have been recently published in j25| under the more 
restrictive assumption on BMO smallness with respect to all independent variables. 

Our approach in this paper is based on the method of approximation which 
was developed by Caffarelli and Peral in [9] in the context of maximal function 
technique, and later adopted for operators with discontinuous coefficients and non- 
smooth domains in [Tll^. We will use maximal functions, the Vitali type covering 
lemma and scaling arguments in L^-estimate for the gradient of the weak solu- 
tion to (jl.ip . in order to derive suitable decay estimates for the level sets of the 
Hardy-Littlewood maximal function for \Du\'^ to increasing levels. After that, the 
standard procedure of weighted integration over the level sets gives the desired es- 
timate ()1.8p . A key point of this approach is comparison with the solutions to a 
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limiting problem, obtained from (jl.ip by averaging in the variables with respect to 
which the coefficients have small BMO norms. The new coefficients depend on the 
remaining one variable, and even if these are only measurable, it turns out that 
the limiting problem supports Vl^^'°°-regularity. 

We would like to point out that it is also possible to use a very influential method 
by Acerbi and Mingione in [TJ El [27]. 

The paper is organized as follows. In Section [2] we introduce the regularity 
assumptions on the coefficients a*^ and the boundary of the domain Q in order 
to state the main result Theorem 12.41 In Section [3] we establish local Lipschitz 
regularity for the solutions to the limiting problem mentioned above. Section|4]deals 
with deriving of L^-gradient estimates from an appropriate perturbation theory in 
analysis. In Section [5] we obtain the optimal gradient estimate (jl.Sp in weighted 
Lebesgue spaces regarding the Dirichlet problem For a particular choice 

of the weight w and as an outgrowth of our main result, we obtain in Section [S] 
gradient estimates in the framework of the Morrey spaces which imply global Holder 
continuity of the weak solutions to generalizing this way the celebrated results 

by De Giorgi [13] and Morrey [28]. 

2. Main result 

We start with the following notations: 

(1) The open ball in R"~^ with center y' = (yi, • • • , Vn-i) and radius r > is 
denoted by 

Bliy') = {x' = {x,,--- , G M""! : \x' - y'\ < r}. 

(2) The cylinder in R"~^ x M with center y — {y',yn) and size r > in the 
x„-axis is denoted by 

Criy) = Bl{y') X (y„ - r,?/„ + r). 

If the center is the origin, we do not specify it and write just Cr for the 
sake of simplicity. 

(3) For each fixed a;„ € M and for each bounded subset U' of R"~^, the integral 
average of a function g(-,Xn) with respect to x'-variables in U' is denoted 

by 

Vu'i^n) = f g{x\xn)dx'^— g{x',Xn)dx', 
Ju' l"^ I Ju' 

and \U'\ stands for the (n — l)-dimensional Lebesgue measure of U' . 

We now state the main assumptions on the data of problem (|l.ip regarding the 
coefficients matrix A{x) and the domain f2. 

Definition 2.1. We say that {A,n) is {S, R) -vanishing of codimension 1 if for 
every point y G and for every number r € (0, R] such that 

dist{y,dn) = min dist{y,x) > V2 r, 

there exists a coordinate system depending on y and r, whose variables we still 
denote by x — (x',x„), such that in this new coordinate system y is the origin and 
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while, for every point y G il and for every number r G (0, R] with 
dist{y,dfl) — min dist{y,x) — dist{y,xo) < \pir 

for some xq G dfl, there exists a coordinate system depending on y and r, whose 
variables we still denote by x — (x' , Xn), such that in this new coordinate system xq 
is the origin, 

(2.2) Csrf^ix: Xn>3rS} C Cs^ n C Cs^nlx: Xn > SrS} 
and 

(2.3) / \A{x',Xn)-AB'^ixn)\''dx<S^. 

Remark 2.2. 1. By a scaling invariance property (see Lemma [575] below) . one can 
take for simplicity i? = 1 or any other constants bigger than 1. On the other hand, 
J is a small positive constant, being invariant under such a scaling. 

2. If {A,il) is ((5, i?)-vanishing of codimension 1, then for each point and for 
each sufficiently small scale there is a coordinate system so that the coefficients 
have small bounded mean oscillation (BMO) in x'-directions with no regularity 
conditions required with respect to the a;„-variable. Regarding the boundary of the 
domain, it is sufficiently flat in the Reifenberg sense in this new coordinate system. 
In other words, the codimension 1 {6, i?)-vanishing property of {A, f2) is a general 
enough condition which is surely satisfied in the particular cases of continuous or 
VMO coefficients a^^{x) and or Lipschitz continuous boundary dil with small 
Lipschitz constant (cf. [331 [M])- It is clear that assumptions (12. ip and (|2.3p allow 
quite arbitrary discontinuities of a*^ (x) in one direction, whereas the discontinuities 
with respect to the remaining variables are controlled in terms of small BMO (think, 
for example, for small multipliers of the Heaviside step function). Moreover, the 
Reifenberg flatness (|2.2p extends the '^-regularity of (jl.ip to the case of domains 
with rough boundaries of fractal nature (cf. |41)). 

3. The Reifenberg flatness condition (|2.2|) implies that the boundary dil satisfies 
the so-called (A)-condition (see [531 [TU]). Namely, setting Br{xo) for the ball of 
radius r and centered at xq, there exists a positive constant Kq(S) such that the 
Lebesgue measure of Br{xo) n is comparable to that of Br{xo) : 

(2.4) Kn{S)\Br{xo)\<\Br{xo)nn\<{l-Kn{S))\Br{xo)\ for each xq G 917. 

4. The numbers V2r and 3r above are selected for our purpose. The reason for 
this selection is that we need to take the size of a cylinder Cr{y) large enough to 
contain its rotations in any directions. 

Before stating our main result, let us recall the definition of the Muckenhoupt 
classes Ag, 1 < s < oo, and the respective weighted Lebesgue spaces L^{il). A 
positive locally integrable function w on R", w G Lj^^CR."^), is called to be a weight. 
Then, given s G (1,00), this weight belongs to the Muckenhoupt class As if 

[w]s = sup sup I + w{x) dx \ i + w{x)^^ dx | < 00, 
y={y'.y^)eti" r>o \JcAy) J yc'Ay) J 

where -f is the integral average and the supremum is taken over all cylinders 

Criy) = {x ^ ix',Xn) G M": \x' - y'\ < r, \xn - yn\ < r}. 
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We note that the As classes are nested, that is, As-^ C Ag^ if 1 < si < S2 < cxd. To 
give an example, consider the function 

u,„(a;) = l^l", xeW\ 

Then Wa G A^ if and only if —n < a < n(s — 1). Thus, Wa is a typical weight which 
can be considered in the present paper. 

For each measurable set E C M" and a weight w, we set 



j{E) = / w{x) 
Je 



dx. 



In what follows, we will use the following important properties of the Ag weights. 

Lemma 2.3. ([40 ) Let w £ A^ for some 1 < s < oo, and let Cr{y) he the cylinder 
Cr{y) centered at y = (y', t/„) G il and of size r > 0. Then we have 





1 

71 



\Criy) n n\ 
\Criy)\ 



< 



w{Criy)nn) 



< 



71 



\Criy)nn\ 



\Criy)\ 

where 71 and (3 > are constants depending only on [w]s and n. 

We next introduce the weighted Lebesgue spaces under consideration in this 
paper. Given a weight w £ As, 1 < s < 00, the weighted Lebesgue space Lf^{^) is 
the set of all measurable functions h: -> M satisfying 



\h\\ 



\h{x)\^ w{x) dx \ < 00 



We are in a position now to state the main result of the paper. 

Theorem 2.4. Given a number p e (2, 00) and a weight w € Ar, there ex- 
ist a small positive constant S — S{i>, L,n,p, [w]e,^) and a positive constant 
c{i>, L,n,p,[w]p,fl) such that if (A, il) is (S, R) -vanishing of codimension 1 and 

\F\ e Lw{fl), then the unique weak solution u € HQ{fl) of satisfies jDup G 

Lw{^) with the estimate 



(2.5) 



Du\Pw{x) dx<c \F\Pw{x) dx. 
n Jn 



Let us point out that the approach employed in the paper is applicable both to 
equations including lower-order terms and to elliptic systems. In order to fix the 
ideas and to avoid unessential technicalities, we limit ourselves to the equations of 
principal type as the one considered in 



3. LiPSCHITZ REGULARITY FOR A LIMITING PROBLEM 

In this section we will prove local Lipschitz regularity for a limiting problem 
which guarantees a fundamental step to derive (j2.5p . For, consider the following 
elliptic equation with coefficients depending on one spatial variable, say x„ : 

(3.1) A {a''ixn)D,v) - in ^2 = {(x',x„): \x'\ < 2, |x„| < 2}. 

Of course, the coefficients a'-' are assumed to satisfy the basic structure conditions 
(|1.2p and (|1.3|) . but these are allowed to be only measurable. Indeed, the solutions 
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under consideration are defined in the weak sense as usual. Namely, we say that 
V g H^{C2) is a weak solution of p.ip if 

/ a''WjvDi(p dx = 

JC2 

for all ifi e H^{C2). 

Throughout this section we denote by c a positive constant that can be computed 
in terms of known quantities such as v, L and n. The standard local i^-estimate 
for (PH is 

\\Dv\\l2^Ci) < cMl^C2)- 
Since the coefficients in depend only on a;„ and the equation is linear, one 
is allowed to differentiate p.ip up to any order with respect to x' — {xi, • • • , a;„_i) 
variables. This observation gives the following lemma. 

Lemma 3.1. If v is a weak solution of (j3.ip . then so are D'v = D^'V with the 
estimate 

\\DD'v\\l2(^Ci) < c\\v\\l2(C2)- 
Moreover, we have the following improving-of-regularity result. 

Lemma 3.2. If v is a weak solution of (j3.ip . then Dv belongs to (Ci) with the 
estimate 

\\Dv\\l2*(^Ci) < c\Ml^C2): 
where 2* is the Sobolev conjugate of 2, that is, 2* = if n > 2 while 2* is an 
arbitrary large number if n — 2. 

Proof. Let u be a weak solution of ((3?T|) . Then by Lemma l3Tl DD'v G L^(C2) 
with the estimate 

WDD'vWl-^^Ci) < c\Ml^C2), 
whence, the Gagliardo-Nirenberg-Sobolev inequality implies that 

(3.2) \\D'v\\L-^c,)<c\\v\\mc2)- 

We next want to show that DnV belongs to (Ci). To do this, remember first 
of all that the coefficients a*-' depend only on Xn and this rewrites the equation 
p.ip in the form 

(3.3) £>„ I a"-" DnV + ^ a'^WjV ] = D^{a'WJv) 

= -'^a'WijV in C2 

in weak sense. Defining 

(3.4) <tJv) = a^DnV + ^ a"JDj W, 

it follows from pT^ . Lemma [O and that 

A(fT„(«)) = a""i?„,z; + ^ a'y'A.^ e L'iCi) 
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for each i ^ n with the estimate 

(3-5) I|A(o-«(w))||l2(Ci) < c\\v\\l2(c^). 

On the other hand, (|1.2p . Lemma IXTl and (I3.3P yield 

and 

(3-6) \\Dn{(yn{v))\\LHCi) < c\\v\\l^C2)- 

It follows from (|33|) and that 
with the estimate 

\Wn{v)\\m{Ci) < c|1w|1l2(C2), 

and applying once again the Gagliardo-Nirenberg-Sobolev inequality we get 

(3.7) an{v) e L^' {Ci) 
with 

lkn(f)||L2.(c^) < c\\v\\h^C2) < c||t'||L2(C2)- 

On the other hand, we have 

by dSll) and (Hill), jS^) and give 

We use the basic structure conditions (11.21) and (11.31) to obtain that D^v G 
L2*(Ci) with 

(3.8) \\DnV\\L2'(Ci) < cM\l^C2)- 

The claim follows from p.2p and p.8p and this completes the proof. □ 

We will prove now interior Vl^^'°°-regularity for the solution of p. II) by employing 
an iteration argument. 

Lemma 3.3. If v is a weak solution of p.ip . then v £ W^'°^ {Ci) and there exists 
a constant c, independent of v, such that 

(3.9) \\Dv\\L^^Ct)<c\\v\\L2(C2)- 

Proof. Let w be a weak solution of p.ip . Then it follows from Lemma [01 that D'v 
are also weak solutions of p.ip . According to Lemma [3.21 DD'v G (Ci) with 
the estimate 

The same arguments as in the proof of Lemma [3?2l applied to 2* instead of 2, give 
Dv G H^i'2*(Ci) and 

\\Dv\\w^,2' iC) < c||w|li2(c,). 
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To proceed further, we consider first the case of an odd dimension n. Repeating 



k times the above procedure gives 2* ' ' ' * — ^" , . Taking k so large that 



n-2k ■ 

k til 



k > 2—1 ensures 2* ' ' ' * > n and we apply the Morrey inequality to find that 



for some < 71 < 1. This implies the estimate p.9p when n is odd and we are 
done. 

Alternatively, if n is even, it may happen that 



n-2k 



after k iterations. However, using the quantity ^_2(^_^\^ — £ instead of for 
sufficiently small £ > 0, we will have 

2n Y f 2n \** 

— £ <n, \ — — ■ - £ > 71, 



n-2{k~l) J ' \n-2{k-l) 
and the proof completes as above. □ 

Our next step consists in proving W^'°°-regularity up to the boundary for the 
weak solutions to the Dirichlet problem for the equation p.ip . To do this, we define 

(3.10) C+ = C2 n {xn > 0} and T2 = = C2 n {x^ = 0} 

and suppose that v € H^{C2) is a weak solution of 

{D,{a^^{x^)D,v)^ inC2+, 
\ t; = on Tz- 

This means that for all ip G i7g(C^), 

a^^ DjVip^ dx = 

and the zero extension of v belongs ro H^(C2)- 

We are ready now to give a natural extension of the interior W^^'°°-regularity in 
Lemma 13.31 up to the flat boundary using a proper reflection argument. 

Lemma 3.4. If v is a weak solution of p. lip , then v S W^'°°{C^) and it satisfies 
the estimate 

||£'«|lioo(c+) < c\\Dv\\^2^c}^ 
with a constant independent of v. 

Proof. Define v in C2 by 

v{x',Xn) if Xn > 0, 



v{x,Xn) 



-v{x' ,—Xn) if Xn <0 



and extend a*^ (a;„) from {x„ > 0} to {xn < 0} by even or odd reflection, depending 
on the indices i and j, in a way that the extended w is a weak solution of (|3.ip . We 
then apply Lemma [3731 to find that v supports the interior VF^'°°-regularity in C2. 
The claim follows from the restriction of v from C2 to C2 ■ Q 
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4. L^-ESTIMATES FROM AN ARGUMENT BY PERTURBATION 

In this section we study gradient estimates of the weak sohition of by 
comparison with solutions to the hmiting problems p.ip and p. lip . The idea is to 
use higher integr ability results, see Lemma 14.11 and Lemma 14.41 These regularity 
results follow from the fine works in the recent papers ^ i30» iSU ES] by Bogelein 
and Parviainen where the authors investigated self improved regularity near the 
boundary with a very mild condition, so-called, the capacity density condition. 
Needless to say, our Reifenberg flat domain satisfies this capacity density condition. 

We will employ here the reverse Holder inequality, as used in [U [2l [6l |8] which 
gives a better regularity of solutions. This reverse Holder inequality can compensate 
the lack of compactness of weak solutions which was previously used in 7 . Thus 
the present approach can be applied to a more general setting when the coefficients 
belongs to with respect to one of the spatial variables where any compactness 
fails. 

We start with the interior case. To do this, let us suppose that F e L^(ri,M") 
and u € H^{n) satisfies 

/ a}-' DjU^if dx — pDiip dx 
Jo. Jn 

for aU if e H^i^l), that is, It is a weak solution of 

(4.1) A {a'Wju) = D,f in fl. 

By a proper translation, scaling and normalization (see Lemma 15.51 below) we 
further assume that 



(4.2) C5 = {{x',xn): < 5,|x„| < 5} C r!, 

(4.3) / \Du\^ dx < 1 



and 



(4.4) / {\A{x',x^)-Ab'S^^)\^ + \F\') 



dx < 6^ 

4 

for some small S > 0. 

Then we consider a local homogeneous boundary problem correspondmg to (|4.ip . 



(4.5) 

and the limiting problem 
(4.6) 



Di{a^iDju) =0 in Cg, 
u = u on 9C3, 



Di a^:>B'.(xn)D,v) =0 in C- 



-■2, 



[ i; = M on dC2- 

We will use the following higher regularity result for the weak solution to (j4.5 



Lemma 4.1. Ifu is the weak solution of (|4.5I) and (|4.2p - (|4.4l) hold true, then there 
exists a constant a = cr(n, i^, L) > such that 

\\Du\\l2+^^C2) < c- 
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Proof. It follows from well-known higher integrability results for (|4.5|) that there 
exists a constant a — (j{n, v, L) such that 



2 + 



(4.7) / \Du\^+'' dx<c[i \Du\^ dx 

JC2 \JC3 

Then the maximality propertj0 and (j4.3p imply that 

(4.8) / \Du\^ dx < c-f iDup dx<c 

and the claim follows from (14.71) and (14. 8p . □ 



Lemma 4.2. If v is the weak solution of (j4.6[) and (|4.2p - (14.41) hold true, then v is 
locally Lipschitz continuous with the uniform bound 

\\Dv\\l^^c,) < c. 

Proof. By Lemma [3.3[ the maximality property for (|4.5p and (|4.6p . and (|4.3p we 
have 

\\Dv\\\^(Ci)^^T \Dv\^dx<c-i \Du\^dx<c+ \Du\'^ dx < c. 

□ 

Lemma 4.3. Given a weak solution u of (j4.ip satisfying (|4.2p and (|4.3p . Zef w 
&e i/ie weafc solution of (|4.6p . T/ien /or < £ < 1 fixed, there exists a small 
5 = 5{e) > SMc/j i/ia< i/ (|4.4p /loWs /or suc/i (5, then we have 

\D{u - v)\'^ dx < e^. 



Proof. Let u be the weak solution of (|4.5I) . Then it follows from a direct computa- 
tion that u — u E HqICs) is the weak solution of 

(4.9) D,{a'Wjiu-u)) = D,r in C3. 
We also see that u — v E Hq{C2) is the weak solution of 

A (^i3^(a;„) A = (^[a'^ {x' ,Xn) - '^Bi{Xn)]DjU^ in C2. 

It follows from the standard L^-estiniate for (|4.9p and the smallness condition (|4.4p 
that 

(4.10) / \D{u-u)\'^ dx<ci \F\^dx<c5'^. 



^Hereafter, by "maximality property" we mean the standard L^-estimate for the weak solutions 
to Dirichlet problem for second-order, divergence form linear elliptic equations. 
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Furthermore, in view of Holder's inequality, Lemma (|1.2p and (|4.4|) we have 
|L»({t-w)p dx= f \A-Ab'Jx„)\^ \Duf dx 



i 



<c[4' U-Ar/(x„)| ' dx\ [j \Du\^+'' dx 







4 + 2<r \ 






(x„) " dx\ 










As', 


(a;,i) dx)j 









< 

< c S"^ 

for some di = (Ti(n, L, v) > 0. In the hght of (|4.10p . we thus deduce 

\D{u - v)\^ dx < c{S^ + 6"')= e^, 

by taking (5 > so smaU, in order to get the last equahty. This completes the 
proof. □ 



We next extend the interior estimates from Lemma 14.31 to obtain boundary gra- 
dient estimates in Reifenberg domains. Recalhng the notations from Section [2] and 
p.lOp . we add here some more geometric notations: 

(1) Qr = flnCr, C+ = Crn{x ^ {x',Xn): X„ > 0}, 

(2) d.^nr ^Onn a, Tr = Crn{x= (x', Xn) ■■ Xn = 0}. 

Motivated from our geometric assumption in Definition 12. 1[ we further assume 
that there exists a small S > such that 

(4.11) Cflr, cnr\{x^ {x,Xn). Xn>-lOS}, 
and 

(4.12) £ (^\A{x',Xn)-AB',{Xn)f + \Ff)dx<S^ 

holds true for such small S. Then let us suppose that u is a weak solution of 



(4.13) 

which means that for all (p E iJQ(ri5) 



D,{a'W,u) = D,r in f^g, 

u = on dyjVt^, 



I a^-' DjuDiip dx = / f^Diipdx 

•/05 J Us 



and the zero extension of u belongs to Hq{C5). By means of suitable scaling and 
normalization, we may also assume that 

(4.14) / \Du\^ dx < 1. 

Then consider the following homogeneous problem 

I w = u on 9515. 
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We then look at the foUowing hmiting problems 

i D, (d:^g,^{xn)Djhj = in f24, 



(4.16) 
and 
(4.17) 



h = w on dw^A, 

Di (^g,^{xn)Djvj = in , 
V ^ on T4. 



Lemma 4.4. Let w be the weak solution of (|4.15p satisfying (|4.1ip . (|4.12p and 
(|4l4| . Then 

\\Dw\\L2+.2(^2i) < C, 

for some small universal constant > 0. 

Proof. It follows from the Reifenberg flatness condition (|2.2p (cf. Remark 12.21 3) 
that satisfies the (A)-condition, that is, the Lebesgue measure of ^r{y) is S- 
comparable to \Cr\ : 

■ 1 _ A"l " 

\Cr{y)\, Vyedn, Vr>0. 



|a(y)| > my)] > 



2V2 . 

By using this property, one can check that R" \ satisfies the uniform 2-thickness 
condition from [5l l30if3ni32| and, as consequence, the uniform capacity density con- 
dition. Then, according to the improving-of-integrability result up to the boundary 
(cf. [301 ISIl 132]) we have 

/ \Dw\^+''^ dx<c(i \Dw\^ dx] ' <c (J \Duf dx] ^ <c 

for some positive constant (T2 = cr2{n, v, L), after using the maximality property for 
(|4.15p and the assumption (|4.14p . □ 

Based on weak compactness argument, we compare the weak solution (|4.16p with 
a weak solution of (|4.17l) . 

Lemma 4.5. For any e > 0, there exists a small S = 5{e) > such that if h is the 
weak solution of (I4.16P with 

/ |L>/ip dx < 1, 

and (|4.1ip - (j4.12p hold for such d, then there exists a weak solution v of J^.i?] ) such 
that 

(4.18) / \Dv\'^dx<l and / \h - v]"^ dx < e"^ . 

Jet Jet 

Proof. If not, there would exist Sq > 0, {hk}^^i and {ri|}^^ such that hk is a 
weak solution of 



A {aijg, ixn)Djhk ) = innl, 



hk = w on dw^4 



(4.19) 
with 

(4.20) C5+cr!gcC5n|x„>-y|, 
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Dhkl"^ dx < 1, 



,^{xn)Djvj = in C^, 
V — on T4, 



J- \Dv\^ dx < 1, 
Jet 



\hk - dx > el- 

We now extend by zero from Vl\ to C4 and denote it by hk also. Then it 
follows from Poincare's inequality and (I4.2ip that ||/ife||_f/i(C4) 1^ c. Then there 
exists a subsequence, which we still denote by {/ife}, and ho S H^{C^) such that 

(4.25) hk ho weakly in H^{C^) and hk ho strongly in L'^{C^). 

From (|4.19p . (|4.20p and ()4.25p we see that ho is a weak solution of 

A (ai7B^(a;n)£'i^o) = in C^, 
/iQ = on T4. 

We next observe from ()4.2ip , ()4.25p and weak lower semicontinuity property that 
i \Dho\^ dx < liminf f \Dhk\'^ dx < 1. 

We then reach a contradiction to (|4.24p from (|4.25p - (|4.26p . This completes the 
proof. □ 



(4.26) 



Now we are in a position to obtain an analogue of Lemma 14731 regarding gradient 
estimates up to the boundary for the weak solutions of (jl.ip . 

Lemma 4.6. Let u be a weak solution of (|4.13p . Then for any £ > 0, there is a 
small S — (5(e) > such that if (^777]), j.^. j|j| ) and ^4-14^ hold, then there exists a 
weak solution v of ^■l'7\ l such that 

(4.27) \\Dvo\\L^<^n,)<c 

and 

\D{u^vo)\^ <e^, 



where vo is the zero extension of v from to C4 . 

Proof. Let w be the weak solution of (|4.15p and h the weak solution of (|4.16l) . Then 
we can derive in a similar way as in the proof of Lemma 14.31 that 

(4.28) / \D{u~h)\^ dx<cS^. 
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From the maximality properties for ()4.15p and ()4.16|) . we obtain 

(4.29) / \Dh\^ dx<c \Dw\^ dx<c \Du\^ dx < c. 
jQi Jn^ 

In light of (|4.29l) and Lemma [4.51 we see that there is a weak solution v of (|4.17p 
such that 

(4.30) / IDwp dx<c 

■let 

and 

(4.31) / \h - dx < c^e^ , 



where c* is to be determined small in a universal way. 

We extend v by zero from to C4 and denote it by vq. Then we derive from 
Lemma EH and that 

(4.32) ||D«o|Uc.(03) = ||i?HlL=o(C3+) < c\\Dv\y(cr) ^ ^' 

which is ([427)1 . 

A direct computation shows that vq is a weak solution of 



Di{aijg,^{xn)D.jV^^ = -Dng"^ in Vt^, 
iiQ = on dw^i, 



if Xn > 0, 

^nn B'. (*^n 

)DnV{x' ,Q) if Xn < 0. 



Di(^aijg,^{xn)Dj[h-vo]j = D^g" in fl^, 
h — vq ^ on dwfl4. 



(4.33) 
where 

(4.34) g'- 
Then in light of (I4.16P and (|4.33p , we find that h — vq is a, weak solution of 

(4.35) 

From standard i'^-estimate for (|4.35p . we have 

(4.36) / \D{h^vo)\^ dx<c( f {\h~vo\^ + \g''\^)dx) . 
We estimate the right-hand side of (|4.36p as follows: 

(4.37) / \h-vo\^dx=~f \h - v\'^ dx + [ \h\^ dx 
Jna Jet l"3| JnsXct 

<c,e^ + j^( f \h\^ dx) " \n3\C+ 



1^31 \Jn^ 
<c.,e^+c6^f \Dh\^ dx 

Here in the first line, we have used (|4.1ip and the fact that vo = in \ . In 
the second line, we have used (|4.3ip and Holder's inequality. In the third line, we 
have used (|4.1ip and Sobolev inequality. In the last line, we have used (|4.29p . 
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(4.38) ^ \9''?dx<-^f a;^B'Ax„)D„v{x',0)) 



Using dual), (jilTI) . (1432]) and (jiTM)) . we deduce 

2 

dx < c6. 

But then (I13S1), (|i37|) and imply 

(4.39) / \D{h - vo)\^ dx < c (c^e^ + 6^) . 

Jn-2 ^ ' 

We now combine (|4.28p and (|4.39p . to obtain 

\D{u - t;o)P dx<c{b^ c.e^ + < c (c^e^ + <5*) . 

Finally, we take c* and (5 so small, in order to get the required estimate 

i \D{u-vq)'^ dx<e^. 
Jn2 

This completes the proof. □ 

5. Gradient estimates in weighted Lebesgue spaces 

In this section we will obtain the optimal weighted W^^'P-regularity for the Dirich- 
let problem (jl.ip . The main analytic tool of our approach is the maximal function, 
so let us recall first of all its definition and basic properties, see [Ml |37] . 



Definition 5.1. Given a locally integrable function h € Lj^^lW"), its Hardy-Little- 
wood maximal function Mh is given by 



{Mh){x) = sup f \f{y)\ dy, 

Cr(x)JCr(x) 

where the supremum is taken over all cylinders Cr{x) in R" centered at the point 
X and of size r > 0. If h is defined only in a bounded domain U C M", we define its 
restricted maximal function as 

Muh = M (hxu) , 

where 

I h{x) if X ^ U, 



{hxu){x) 



ifx^U. 



The important properties of the maximal function with respect to weights are 
summarized in the following 

Lemma 5.2. (|:40J Given a weight w ^ for some s e (1,cxd), suppose that 
h e L^(M") C LlM")- Then we have 

^ll'i||L-^(R") < II-M^IIl-,(R") < c|1/i|!lj^(R"), 

where c ~ c{n, s, [w]s) > is a universal constant. In the particular case w{x) — 1, 
it follows 

(5.1) |{a;eR": (Mh) (x) > X}\ < cj J \h{x) dx, (A > 0), 

with a universal constant c — c{n) > 0. 
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We will use also the following standard result from the classical measure theory 
regarding weighted Lebesgue spaces. 

Lemma 5.3. ([40 ) Given a weight w € As for some s G (l,cx)), let h be a non- 
negative function lying in L^jiJJ) C L^{U) for some bounded domain U C M". Let 
9 > and X > 1 be constants. Then 

h e LliU) ^ S^J2 ^'"'^{{^ € U: h{x) > ex''}) < oo 

k>l 

with 

ls<\\h\\l.^^ij^<ciw{U) + S), 

for some universal constant c — c{0, A, s). 

In what follows, our approach is mainly based on the following version of the 
Vitali covering lemma stated in the settings of weighted measurable sets. 

Lemma 5.4. Let fl be a bounded, ((5, 1)-Reifenberg fiat set, and let w be a weight 
from As for some s € (l,oo). Suppose the measurable sets D C E C ft have the 
following properties: 

(5.2) w [Df] Ci(y)) < £w (Ci{y)) for some e £ (0, 1) and for every y ^ fl; 



(5.3) 

Then there holds 



for any cylinder Cr{y) with r > and y £ 51, 
w{D n Criy)) > ew{Cr{y)) implies ^ n Cr(jj) C E. 



w{D) < 72e w{E) 

for some positive constant 72 = 72 ((5, n, s, ['w]s). 

Proof. Fix any y € D and consider a continuous function 

, , w(DnCr(y)) 
w{Cr[y)) 

Then w{Q) — \i'a\r^o+ 'cci{r) — 1. From the assumption (|5.2|) . we see w{\) < e. 
Thus, one can find a number ry = r{y) G (0, 1) such that 

(5.4) ^{^y) — £1 ■uj{r) < £, Vr > ry. 

Since {D n Cr^ iy)}ye d is an open covering of D, there exist a disjoint subcovering 
{Criiyi)}i>i with ri = r{yi) < 1 such that 

(5.5) \J[DnC5rM)]=D. 
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We now compute as follows: 



w{D) 



< 



JO! 

< 



< 



^i>l J i>l 

i>l 

c£ w (Cniyi)) by the doubling property of w, 



i>i 



Lcmina l2.3l 

< 



(ED 

< 



C£ 



E 

2V2 



\CrM)\ 



10 



nnCrAy^)\ 



2^/2 



n;3 



e w I [J[r2 n Cniyi)] ] by the disjoint covering, 



< 



2V2 

1-S 

This completes the proof. 



-1 710 



ew{E) 



□ 



We will employ also the following invariance property under scaling and normal- 
ization which follows by a straightforward computations. 



Lemma 5.5. Under the scaling 

Ap{x) — A{px) and Vtp = I —x: a; e 



(0<p< 1), 



the structure assumptions (jl.2p . (|1.3p and the main assumption in Definition \2.1\ 
are invariant with the same constants v, L, and the dilated scale —. Furthermore, 
along with the normalization 

ux,p{x) = and Fx,p{x) = — (A > 1), 

ux^p{x) is a weak solution of 

I A {a'^{x)D,ux.p{x)) - D,flp{x) m ^p, 
1 ux.p{x) =0 on dflp. 

Based on the maximal function and scaling argument and making use of Lemma f4.3l 
and Lcmma l4.61 we have the following result. 

p 

Lemma 5.6. Let |Fp e Li,{n) e L^{rt) with w e Ap, 2 < p < 00, and let 
u G Hq{VL) be the weak solution of (jl.ip . Then, there exists a universal constant 
A2 > 1 such that for each < £ < 1 one can select a small 5 = die) > such that 
if {A,fl) is [S, R) -vanishing of codimension 1 and Cr{y) satisfies 

(5.6) w{{x€n: Mi\Du\'^) > Xl} n Criy)) > s w (Criy)) 



WEIGHTED L" ESTIMATES FOR ELLIPTIC EQUATIONS 



19 



for such a small 6, then we have 

(5.7) Cr{y)nn c {xen-. m{\Du\^) > i}u {x en-. m{\f\^)>s^}. 

Proof. We argue by contradiction. Assume that Cr{y) satisfies (|5.6I) but the claim 
(|5.7p is false. Then there exists a point yi e ^r{y) = Cr{y) H ^ such that for each 
p > one has 



(5^8) 



\Du\^ dx < 1 and i \F\'^ dx < 5^. 



We first investigate the interior case when C'g^^(y) C 51. According to Defini- 
tion 12.11 there exists a new coordinate system, modulo reorientation of the axes 
and translation, depending on y and r, whose variables we denote by z, such that 
in this new coordinate system y = 0, j/i — zi^ 



(5.9) 
and 



/ A{z' 



'C'4%/2r 

From ()5.8p and ()5.9p we see that 

'5 



^n) - ^B' ^ (Zn) 

4 V2r 



L 



IDur dz < 



and 



iFr dz < 



Applying Lemma [5.51 with p = \f2r and A = yd)", the last three inequalities 
imply that we are under the hypotheses of Lemma 14.31 which gives, after scaling 
back, that for the weak solution v of 



(5.10) 

one has 
(5.11) 

and 
(5.12) 



= in C2^^, 



on dC. 



2\f2r 



\D(u~v)f dz < 



V2t 



||i?«||L^»(C^„)<Ao 



for some positive constant Aq — Xo{iy, L,n). 
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We have 

{zen-. M{\Du\^) > 8X1} n Cr 



[zeQ: Mc,^^i\Diu-v)f)>Xl}nC^^ 
+ \{zen: Mc,^^i\Dv\')>Xl}nC^^ 

^ {zen: Mc,^^mu~v)\^)>xl}nc^,^ 

\D{u-v)f dz 



15T1 

< c 

tsTTTl 

< coe^iai, 

for some positive constant cq = cq(i>, L, n). In other words, 

\{x£n: Mi\Du\^) > 8X1} n a(y)| < coe^\Criy)\, 
whence Lemma 12.31 vields 

(5.13) w{{xeCr{y): M{\Du\^) > 8X1}) < (cos^f w {Cr{y)) 

with a positive constant ci = ci ^z^, L, n, ■ 

We next consider the boundary case when Cr^^^{y) ^ fl. Now, for the sake of 
simphcity, we denote c to mean a universal constant c = c(z/, L, n) that is indepen- 
dent of S. We may also suppose that there is a boundary point yo G Sfi fl '^^^^(j/). 
According to Definition 12.11 there exists a new coordinate system, modulo reorien- 
tation of the spatial axes and translation, depending on yo and r, whose variables 
we denote by z, such that in this new coordinate system the origin is yo + SoHq for 
some small Sq > and some inward unit normal nf) to dfl ai yQ, y — zq, yi — zi, 

(5-14) ^^2^472. ^ ^24v^,. C {z e C,,^^: z„ > -48V2rS} , 

and 

(5.15) / A{z',Zn 

Then it follows from (|5J4)) and ((5?T5l) that 

(5.16) / \Dufdz<2(^y-f 



24v'2r 



and 
(5.17) 



/ 
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|Z?u|" < 2 ( — 

" 24 



< 2 ( ^ 1 52 



Applying Lemma [531 with p = 6\/2r and A = -\/2 (||)", ([g?Ti)) - (ICT7)) show that 
we are under the hypotheses of Lemma l4.6l which gives after backscaling that there 
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exists a weak solution v of 



(5.18) 

such that 
(5.19) 

and 
(5.20) 



V = Q on r,. 



/ \D{u~Vo)\ 



12V2ri 

2 < 



where vn is the zero extension of v from r- to Cn /o . 

12\/2r 12\/2r 

We have 



|{zerJ: XdDup) > 8A?}na(zo) 
JsTTH i.i lsTTsl i (- 

1^ e f^6V2^ 



< 



15381 

< 



J530J 



< 



C, 



6V2r 



{zefi: Xo,,^,(|i^(u-^;o)P)>A?}n 



C, 



6V2r 



12\/2t 



< C |f^l2v/2rl < C e\Cr[,ZQ)\. 

This way 

I {x e 1^: XdZJup) > BA^} na(2/)| < c£|a(2/)| 
and therefore Lemma 12.31 imphes 

(5.21) w{{x<^^: X(|L»up) > 8A?} na(y)) < C2e'^u>(a(2/)) 

with a constant C2 = C2 (^v^LjU, [w]e^ > 0. 

Finally, we combine (15 131) and (j5.2ip and set A2 = max{Ao, Ai} in order to reach 
a contradiction with (15. 6p since e is arbitrary given. This completes the proof. □ 



Fix now e and take S and A2 as given in Lemma 15.61 We use the Vitali covering 
lemma (Lemma 15. 4p in order to obtain the power decay of 

w{{{x,t) e ft: M{\Duf) > Xl}) . 

Lemma 5.7. Let |Fp e Ll,{VL) G L^{^) with w G Ar, 2 < p < 00, and let 
u G H^[Vl) he the weak solution of p.ip . Suppose (A, il) is {5, R) -vanishing of 
codimension 1 and set £\ — 7ie. Then we have 

(5.22) w{{x^n: M{\Du\'^)> \f}) <e\w{{x£Vl: M{\Du\^) > \}) 

k 

+ Y.e\w({xen: X(|F|2) >,52^2(fc-»)|^ (/fc = 1, 2, • • • )• 
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Proof. We set 

D^{xe n-. M{\Du\^) > Xl} 

and 

E^{xen: M{\Du\^) > l}U{x en-. M{\F\^)>5^}. 



Clearly, D d E d Q and w (D H Ci{y)) < ew {Ci{y)) by Lemma [2?3l and a proper 
normalization of the problem as one can take A large enough in Lemma 15.51 

The second condition of Lemma [5.4l is a direct consequence of Lemma [5.6l Applying 
Lemma 15. 4i we get the claim in the case k — 1 and iteration of the foregoing 
arguments leads to (|5.22p for each fc > 1. □ 



We are now ready to give a complete proof of our main result, Theorem l2.4l For, 
we will employ Lemma [531 with s — ^, h — A^(|£)up), \ — X2, — 1. We have 



Y^xT-'w{{M{\Du\')>Xf}) 
fc=i 

< ^A2'''^^e^w({.Te ^1: A^d^Mp) > 1}) 

k=l 

00 k 

+ E ^^^^ E {i-^^-- -^(i^n > sX'"-''}) 

00 

fc=i 

00 / 00 

+ E (^2^1)' E^^^'"''^- ({- e f^: MilF^ > S'xf'-'^}) 

i—1 \k—i / 

Lemma l5.3l / ^ p \ ■s^^ i. 

< c(w{n) + \\M{\F\')\\\ E(^2£i) 



k 



Lcmma l5.2l / ^ ^ \ ^ — ^ t. 

for some universal constant c = c (^S, v., L, n, [wjz^ > 0. 

We then first select e > so small to have X^ei < 1. In view of Lemma [5761 one 
can find a small constant 6 = 6 (^^L^n, [w\e^ such that 

00 / 
Y^xf-^^w {{M{\Du\') > Xf}) < c (w{n) + lllFpp, 



fe=i 



lS. (n) 



holds true for now fixed small (5 > for all {A, ft) which are {6, i?)-vanishing of 
codimension 1. Therefore, it follows from Lemmae 15.21 and 15.31 that 



\Du\Pw{x) dx<c\^J \F\Pw{x) dx + w{n) j , 

which implies the desired estimate (j2.5p through the Banach inverse mapping the- 
orem. This completes the proof of Theorem 12.41 □ 
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6. MORREY REGULARITY OF THE WEAK SOLUTION 

We will apply now our main result Theorem 12.41 in obtaining gradient estimates 
in Morrey spaces for the weak solutions to (ll.ip . 

Let us start, first of all, with recalling the definition of the Morrey spaces. Given 
a bounded domain Q C M", p e (1, oo) and A e (0, n), a function / g LP{Q) is said 
to belong to the Morrey space LP '^{^1) if 



ll/IUp.^(n) = sup {-J \fix)\^ dx\ < CX3, 

a;oeO,re(0,diamr2) y JB^(xo)nSl J 

where i?r(xo) is a ball centered at xq and of radius r > 0. The above quantity 
defines a norm under which LP''^{Q) becomes a Banach space. The limit cases 
A = and A = n give rise to LP{V,) and L°°{V,), respectively. 

The following result extends the VF^'^(J7)-regularity theory of to the set- 
tings of Morrey spaces. 

Theorem 6.1. Assume (|1.2p and (|1.3p . Given p S (2,oo) and A S (0,n), there 
exists a small positive constant 6 and a positive constant c, depending onn, L, v, p, X 
and f2, such that if {A, J7) is (6, R)-vanishing of codimension 1 and F G LP'^{^, R"), 
then the unique weak solution u e HqITI) of the Dirichlet problem satisfies 
Du e LP'^(fi,M") with the estimate 

(6.1) ||Du||ip,A(n,E") < c||i^||ip,A(o_R„). 

Proof. Extend F as zero outside O and fix arbitrary xo & V, and r > 0. Set XBr{xa) 
for the characteristic function of the ball Br{xo) and MxBr(xo){x) fo^' its Hardy- 
Littlewood maximal function. 

It follows from Proposition 2 in |T2] that if cr G (0,1) then {MxBr{xo){^)) 
belongs to the Muckenhoupt class Ai, that is, 

M(^{MxBAx„)ix)y) < c {MxBAxo)i^)y for a.a. x e M", 

whence {M.XBr{xo){^))'^ ^ for each p e (2,cx)) with 

[(A^XB.(a:o)(a;))'']p =c{n,p,a). 



Choosing cr e , it follows from Theorem 12.41 that there exist constants S 

and c, depending on n, L, i^, p, A and fl, such that if {A, O) is ((5, i?)-vanishing of 
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codimension 1, then 



(6.2) / \Du{x)\P dx = \Du{x)\P {xBA.o)i^)y dx 

JBr(xo)nn Jn 



< 



< 



\Du{x)\P {MxBA.o)i=^)y dx 



c / \F{x)\P {MxBAx„)ix)y dx 



= c J jFix)\P {MxBAxo){x)y dx 
I 



B2r{X0) 



\F{x)\P {MxBAxo)ix)y dx 



Ioir,xo) 



\ 



E/ \FixW {MxBAxo){x)y dx 



Ikir,xo) 



J 



after the dyadic decomposition M" = i?2r(2;o)U (UfcLi B2k+ir{xo) \ B2kj.{xQ) 
We have MxBr(xo){x) < 1 a.e. K" whence 



(6.3) 



Io{r,xo)< / \F{x)\P dx<c{ny\\F\\l 

JB2AX0) 



To estimate Ik{r,xo), we note that for each x € B2k+ij,{xa) \ B2kr{xo) and each 
p > we have 



(6.4) 



Moreover, having in mind x G B2k+ir{xo) \ B2kj,{xo), the term on the left-hand side 
above is positive only for values of p greater than — r. This way, the obvious 
inequality 2*^' - 1 > 2''-'^ Vfc > 1 reduces dO]) to 



1 



I r" 1 

XBr(xo) \ {x) dx < 2„(A;-l)^n ^ 2"('=-l) 



and taking supremum with respect to p > we obtain 

1 



(XXS.(:.o)(2;)) < 
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Therefore, 

(6.5) i,(^r,xo)<—^f \F{x)\Pdx 

< I \F{^r dx 

1 f \F(x)\P dx 

^ oA+crri/oA— crnNA; An plip 

- ^ ; ^ 11^ llLp.^(n,K")- 

Substitution of (gSl) and into (H^l) yields 




B,^{xo)nn 

with a convergent series thanks to the choice a £ l) . Dividing the both sides 
of the last inequality by and taking the supremum with respect to xq G ft and 
r > gives the desired estimate (|6.1I) . □ 

The gradient estimate (I6.ip and the known properties of functions with Morrey 
regular gradient (see Lemmae 3. Ill and 3. IV in [101) imply immediately better 
integrability and Holder continuity of the weak solution to (jl.ip for appropriate 
values of p and A. Namely, 

Corollary 6.2. Under the assumptions of Theorem \6.1\ let u £ H}f{Q.) be the weak 
solution of (jl.ip . Then 

(1) u e L^'^(n) c LP^^+P{Q.) ifp + \<n] 

(2) u E LP (f2) for any p' < oo and any A' < n, if p + X ~ n; 

(3) u e C^'^-'^m ifp + \>n. 

It is worth noting that the global Holder continuity with some exponent for 
the weak solutions to divergence form elliptic equations with only measurable co- 
efficients is the essence of the celebrated De Giorgi regularity result ([l^) when 
F e LP{Q.,W) with p> n, and that of Morrey [28] when F e LP^^{n,W') with 
p + X > n, both holding in domains with Holder continuous boundaries. Apart from 
the fact that we are dealing with Reifenberg flat domains, in our more restricted 
situation (coefficients which are measurable in one variable and small BMO in the 
remaining ones) we provide an explicit expression for the Holder exponent of the 
weak solution. 

The results from this section will be applied in a forthcoming paper to the 
study of Morrey regularity of weak solutions to quasilinear divergence form elliptic 
equations with controlled growths of the nonlinearities. 
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